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Abstract 

In this paper we show how the perturbative procedure known as stochastic limit may be 
useful in the analysis of the Open BCS model discussed by Buffet and Martin as a spin 
system interacting with a fermionic reservoir. In particular we show how the same values 
of the critical temperature and of the order parameters can be found with a significantly 
simpler approach. 
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I Introduction 



In this paper we analyze the Open BCS model as given in [1, 2] using the techniques of the 
stochastic limit approach (SLA), which is described in many details in the monograph [3]. 
Instead of considering a fermionic reservoir, as the authors do in [1, 2] (following the original 
suggestion contained in [4] which allow to avoid dealing with unbounded operators), we will 
consider here a bosonic thermal bath. This choice is made to try to stay closer to the real 
physical world, where the reservoir is bosonic. This means that some of our formulas are only 
formal, but they can be made rigorous with just a little effort, using, for instance, the same 
framework for unbounded operators developed in [5] and references therein. We will comment 
again on this aspect of our model in the next section. 

The main outcome of this paper is that the same values of the critical temperature and of 
the order parameters can be found using the SLA, in a significantly simpler way, as we will 
show in Section III. This simplification allows us to focus our attention on some aspects of the 
model which could appear not so clearly using the standard technique. This is what has been 
already observed in other physical applications: for instance, in [6], we used the SLA to explore 
in details some relations between different models of matter interacting with the radiation, 
as the Hepp-Lieb and the Alli-Sewell models. Also, in [7], the SLA was used in connection 
with the fractional quantum Hall effect, giving some interesting results. Other applications are 
contained in [3] and references therein. 

The paper is organized as follows: 

in the next section we introduce the model and compute its generator using the SLA together 
with a semiclassical approximation, already introduced in [2] , useful to obtain the free evolution 
of the matter operators; 

in Section III we write the equations of motion for some macroscopic variables of the matter 
and we recover the same results as in [1]; 

our conclusions are contained in Section IV, while the Appendix is devoted to review some facts 
concerning the SLA, useful to keep the paper self contained. 

II The Physicals Model and its stochastic limit 

Our model consists of two main ingredients, the system, which is described by spin variables, 
and the reservoir, which is given in terms of bosonic operators. It is contained in a box of 
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volume V = L 3 , with N lattice sites. We define, following [1, 2] 

N N 
j=l i,j=l 

where the indexes i,j represent the discrete values of the momentum that an electron in a fixed 
volume can have, o1~ creates a Cooper pair with given momentum while a~ annihilates the 
same pair, e is the energy of a single electron and — g < is the interaction close to the Fermi 
surface. As we can see, only the ± component of the spin, that is the x, y components, have 
a mean field interaction, while the z component interacts with a constant external magnetic 
field. The algebra of the Pauli matrices is given by 

[at, aj] = 8 i3 al [of, <r°] = T^af. (2.2) 

We will use the following realization of these matrices: 

If we now define the following operators, 

1 N 

S N = TjYl a i' R n = S^Sx = R N , (2.3) 
JV 1=1 

H^ ys ^ can be simply written as H^ ys ^ = N(eS%—gR N ) and it is easy to check that the following 
commutation rules hold: 

[S%,R N ] = [H^ ys \R N ] = [H^ ys \S%] = 0, 

for any given N > 0. It is also worth noticing that the intensive operators are all bounded 
by 1 in the operator norm, and that the commutators [S^, <Jj] go to zero in norm as when 
N — > oo, for all j, a and (5. 

As we have already mentioned in the Introduction, we consider here a realistic bosonic 
reservoir, so that some of the following formulas must be understood to be formal. However, 
using for instance the same algebraic framework discussed in [5] for some different spin-bosons 
models, or replacing the bosonic operators with their smeared versions, everything can be made 
rigorous. We avoid here this useless complication, since it would make all the treatment much 
more complicated, hiding in this way our main results. 



3 



Our construction of the reservoir follows the same steps given in [2] , but for the commutation 
rules. We introduce here as many bosonic modes a^j as lattice sites are present in V. This 
means that j = 1,2, ...,7V. p is the value of the momentum of the j-th boson which, if we 
impose periodic boundary condition on the wave functions, has necessarily the form p = ^n, 
where n = (ni,n 2 ,n 3 ) with rij G Z. These operators satisfy the following CCR, 

[ap,i, aq,j] = Wp,i, 4j] = °> °ij = S H S P? ( 2 - 4 ) 

and their free dynamics is given by 

H (re s) = J2J2 <,;«:,/<,.,■ (2-5) 

J=lpeAjv 

Here is the set of values which p may take, according to the previous remark: Ajv = {p = 
^-n, n E Z 3 }. It is useful to stress that the energy of the different bosons is clearly independent 
of the lattice site: = = 4n ^l^a?^^ ■ 

The form of the interaction between reservoir and system is assumed to be of the following 
form: 

/4 7) = £«%(/) + ^ c -), (2-6) 

3=1 

where we have introduced %(/) = J2peA N a p,jf(p), f being a given test function which will be 
asked to satisfy some extra conditions, see equation (2.24) below and the related discussion. 
We would like to stress that, in order to keep the notation reasonably simple, we will not use 
the tensor product symbol along this paper whenever the meaning of the symbols is clear. 
The finite volume open system is now described by the following hamiltonian, 

H N = H° N + \Htf, where H° N = H^ ys) + H^ es) (2.7) 

and A is the coupling constant. 

The first step in the SLA is the computation of the free evolution of the interaction hamil- 
tonian: 

N 

H$(t) = e iH >H^e- iH °^ = ^^^e-^^e^'a^e-^ 1 + h.c). (2.8) 

3=1 

The computation of the part of the reservoir is trivial and produces 

-rj-(res), ■TT( res )j. , -x , 

e lH N \-(/)e- l < '* = aj (fe- lte ), 
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where a,j(fe~' lte ) = J2pe\ N a p,jf(p) e ~ ttep - This is an easy consequence of the CCR (2.4). The 
free evolution of the spin operators is more difficult and its expression can be found in [1, 2], for 
instance. Here it is shown how to obtain the time evolution in a semiclassical approximation, 
i.e., when the free time evolution of the intensive operators S% is replaced by its limit (taken 
in the strong topology restricted to a certain family of relevant vectors, [8]). 
The differential equations of motion for the spin variables are 



da + (t) 

da°At) _ _ , - , „ 



, -2iia+(t)+igS+(t)a°(t) 

-2ig(af(t)S^t)-ar(t)S%(t)). 



where we have called, with a little abuse of language which is quite useful to maintain the 

■ tt( s V s ) + ■ tt( s V s ) J 

notation simple, cfyt) = e N °""e~ N ■ in fact, to be more precise, instead of cr°(t), we 
should write aj'//~ ee (t), to stress the fact that e tH N V ^o^e -111 ^ >f only produces the free evolution 
of a", i.e., the evolution without any reservoir, and for N fixed. Moreover, in (2.9) we have 
introduced S%(t) = e^^e" 1 ^* = ± £f =1 e iH ^* ) Vfe""^* = ^Ef=i^ a (0- 

Let us now call S a = T — strong lim^^^ S 1 ^. The proof of the existence of this limit 
(together with all its powers) may be found in [8] and references therein. We can now take the 
sum over j = 1,2, ...,N of (both sides of) the equations in (2.9), divide the result by N, and 
consider the T — strong limjy^oo of the equations obtained in this way. We find that S°(t) = 
and S + (t) = i(2e + gS°(t))S + (t). These equations can be easily solved: S°(t) = S° = (S°Y and 
S+(t) = S + e i( - 2i+9S °^. Of course S~(t) = (S + (0. The system (2.9) gives now, if we replace 
S%(t) with its T — strong limit S a (t), 

dcr + (t) 

daf(t) _ _ ( - IU ' 



(// -2ie<j+(t)+igS + (t)<j°(t) 
,„ ^ - 2ig(a+(t)S-(t)-aT(jk)S+{t)). 



This system is called the semiclassical approximation of (2.9), and it can be explicitly solved 
using, for instance, the Laplace transform technique. The computation is rather long and we 
omit here all the details, which can be found in [1, 2] . Also, since only o~+(t) appear in (2.8), 
together with its hermitian conjugate, we give here only the result we need. We have 

<r+(t) = e l Vo + e %{v+u))t p> + + e^- w V-, (2.11) 

where we have defined the following operators 

7° + 2S + a~) 

+ ^a°-gS+ar) (2.12) 



Pi 


_ s 2 s+ 

~ u> 2 


(2S~ 


-a+ + S° 


p\ 


_ 9S+ 
~ w 2 


(gs- 


u-gS° + 
uj+gS° U j 


P 3 - 


_ 9S+ 
- u 2 


{gs- 


uj+gS° + 
w-gS° a j 
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and the following quantities 

u = g^J{S ) 2 + 4S+S-, u = 2~e + gS°. (2.13) 

Denning further 

Va(j>) = v - £p + OiU, (2.14) 
where a takes the values 0, + and — , the operator H^\t) in (2.8) becomes 

N 

H { N ] (t) = E E {Mfe ltUa ) + h.c) . (2.15) 

j=l a=0,± 

Remark:— It may be worth remarking that we would have obtained exactly this free time 
evolution even for a fermionic reservoir, since CCR and CAR produce the same free time 
evolution for both the annihilation and the creation operators. From this point of view, the 
difference between a fermionic and a bosonic thermal bath appears really only a minor aspect 
of the model. 

The next step in the SLA consists in computing the following quantity 

h(t) = (-£ff dh £ dt 2 u tot (Hj?(^)HjP(^j) , (2.16) 

and its limit for A going to zero. Here the state lOtot is the following product state uj to t = &sys up, 
where uj sys is a state of the system, while up is a state of the reservoir, which we will take to be a 
KMS state corresponding to an inverse temperature f3 — It is convenient here to use the so- 
called canonical representation of thermal states, [3], which is sketched in the Appendix. Then 
we introduce two sets of mutually commuting bosonic operators {c^j }, 7 = a, b, as follows: 

a^ = ^)cf J + ^)c^\ (2.17) 

where 

. 1 g-/3ep 

m(p) = Ma^ay = n(p) = uj^al-a^) = 1 _ e _ 0£f - (2-18) 

The operators satisfy the following commutation rules 

[<$,c$]=Wrf«r, (2-19) 
while all the other commutators are trivial. Furthermore, we introduce the vacuum of the 



operators , $ : 



cg$ = 0, Vp G A N , j = 1,..N, a = a,b. (2.20) 
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Finally, if we define f m (p) = y/m(p)f(p) and f n (p) = y/n(p)f(p), we get 

ajife^) = cf{f m e^) + cf(f n e u ^), (2.21) 

where we have used the usual following notation cf\g) = Y,pe\ N c pj9(p) an d c^p (g) = 
EpeAiv c p] q(p) 2 - Therefore we have 

H ( N ] (t) = J2 E U (cfifme^) + cfifne**-)) + h.c) , (2.22) 

j=lQ=0,± L V 7 J 

and the KMS state top can be represented as the following vector state, as in a GNS-like 
representation: 

u p (X r ) =< $o,X r $ >, (2.23) 

for any observable of the reservoir, X r , since uop is a gaussian state, [3]. This fact, together 
with (2.20) and with the commutation rules (2.19), simplifies the computation of the two point 
function uj tot (^H^\^)H^\^)^j , which, after some algebraic computations, produces 

= E E E {^.(^ t )l/ m (Pl|V^® e -^''^+ 

V 7 j=la,/3=0,±peAjv 

Since we are interested to the limit A — > of I\(t) we need to impose some conditions on the 
test function f(p), [3]. In particular, we will require that the following integral exists finite: 

f dr £ \fr(p)\ 2 e ±iT ^ < oo, (2.24) 
where f r (p) is / m (p) or f n (p) and ^(p*) is given in (2.14). Under this assumption we find that 
7(0 = Km/A(*) = -*E E K S (^aVi a) +^(^Vjri 6) }, (2.25) 



E 

j=l a=0,± 

where the two complex quantities 



rW= /° rfr £ |/m(p)| 2 e-^®, r?= /° rfr £ \f n (p)\ 2 e^ (2.26) 

J—oo -v_ A J— CO — »_ * 



2 It may be worth noticing that both c^\f) and c^\f) are linear in their argument / 
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both exist because of the assumption (2.24). 

To this same result we could also arrive starting with the following stochastic limit hamil- 
tonian 

H?(t) = EE k Ul(t) + <$>\t)) + h.c] , (2.27) 

where the operators c£] (t) are assumed to satisfy the following commutation rule, 

[c%Ht), c$\t')} = 8 jk 5 aP 5„5{t - Or?), for t > if. (2.28) 

We mean that, as it is easily checked, the following quantity 

At) = H) 2 fdti f 1 dmUH^it^H^ih)) 
Jo Jo 

coincides with I(t). Here Q tot = uj sys Q = uj sys < ^ >, where \l/o is the vacuum of the 
operators c^]{t): c^]{t)^Q = for all a,j, 7 and t, [3]. 

Following the SLA, [3], we now use H^ l \t) to compute the generator of the theory. In fact, 
this is the main reason why this operator is introduced in the game. Let X be an observable 
of the system and t r the identity of the reservoir. Its time evolution (after the stochastic limit 
is taken) is j t (X ® l r ) = U}(X £g> l r )U t , where U t is the wave operator satisfying the following 
differential equation d t U t = —iH^ l \t)U t , whose adjoint is d t Ul = ill} (i) . 

Then we find 

d t j t (X <g> l r ) = iU}[H { f(t),X® tr]U t = 

j=la=0,± 1 } 

Next we have to normal order the formula above, i.e. to move to the right all the annihilation 
operators c^]{t) and to the left the creation operators (t). To achieve this result we need 
to compute first the commutator [c*£j (t), U t ], and this can be done easily by means of the time 
consecutive principle, [3], and of the commutation rules (2.28): 

[c$(t),U t ] = -ijy:](t),H^ l \t')]Ut'dt' = -i fyjT^5{t-t'))U t ,dt' = -ip> a ] V^U t . 

(2.29) 

Similarly we get 

[c^(t),U t ] = -ifir^U t , (2.30) 
and, taking the adjoint of these two equations, 

[U!,c$\t)\ = iU}p> a Tg> and [Ulc^{t)\ = iU}p> a W 
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Going back to d t j t (X <g> l r ) we find that 

N f , , s 

d t]t {x®^) = iY, £ { ( tftfpiW + eg \t)un^ a ,x]u t + 

j=l a=0± V ' 

(iul p^rF + c^\t)uf) [fi\x}u t + uM,x] (-iptfrput + u tC $(t)) + 
+f//[pi t ) x] (-ipirw^i + u t c®(tj) } 

which has to be computed on the state Vl tot . Therefore, since the generator L satisfies the 
equality Q tot (d t j t (X <g> l r )) = Q tot (j t (L(X))), we get 

N , s 

l{x) = £ £ xyjv^ + [pi f , x]p> a r® - p> a[pi \ x\rP - ^ V., ^]rL 6) (2.31) 

i=ia=o± L ' 

This expression can be made simpler if the observable X is self-adjoint (X = X^). In this case 
we have 

L(X) = L 1 (X) + L 2 (X), (2.32) 

where 

N N 

M*) = E E {k^R^ + ^.c}, l 2 (x) = E £ {kURrlH^.}. (2.33) 

j=l a=0± j=l a=0± 

This formula will be the starting point for the analysis in the next section. 

Remark:- Before going on, it may be interesting to stress that, when compared with 
the standard perturbative approach for the master equation for open quantum systems, [2] , the 
perturbative approach based on the SLA appears quite friendly. For instance, the so called time 
consecutive principle and the new Hilbert space with ground vector arising after the limit 
A — > is taken, are typical tools of the SLA and they are essential to make many computations 
almost trivial. 

Ill The phase transition 

As discussed in [1, 2], S% and Rn are the relevant variables whose dynamics must be considered 
to analyze the phase structure of the model. These intensive operators are both self-adjoint, so 
that we can use equations (2.32) and (2.33) instead of (2.31). As a matter of fact, in both [1] 
and [2] these equations of motion are considered only as an intermediate step to compute the 
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equation for A N = ^R 1 ^ 2 , which is called the gap operator. We will see in a while that the same 
conclusions as in [1, 2] can be obtained without introducing A at but working directly with Rn 
and S%. 

As a first step we compute L(S%) = Li(S&) + L 2 (S%). We have, using (2.2), (2.3) and 
(2.33) 

i N i N 

M^) = ^EM-?) = ^E E {K^> J M a) + h.c.} 



j=l j=la=0,± 



which can be written as 



Ll(S° N )= E {(fta+^ + fta-^ + ftoO^ + ftail^rW+^C.}, (3.1) 

a=0,± 

where In — j^Y^jLi lj an d the various coefficients {6 a7 } have been introduced here only to 
stress the fact that Li(S N ) is linear in the intensive operators. As we have already mentioned 
before, the limit of the right hand side of the formula exists in the strong topology restricted 
to a certain family T of states, since all the operators S% converge in this topology. Therefore 
also the limit of the left hand side does exist in the same topology. After some non trivial 
algebra we find 

Li(g°) :=F- strong lim Ll {S° N ) = JgyW +KT< B > U + 9 \ 



n^oo 1V Jvy uj 3 { + (w + gS ) 2 ~ (u-gS ) 2 }' 

(3-2) 

where ^RT^ indicates the real part of r±\ 

The computation of L 2 (S°) := T — strong \im N L 2 (S' N ) follows essentially the same steps 
and produces 

MS »> . + ^ } • (3.3) 



so that the final result is 

L(S°) = - 894S ° { ^ +S ~ )2 h(S°,S + S-). 
Here we have introduced, for brevity, the function 



(3.4) 



(3.5) 

and we have made explicit the fact that h depends on S + S~ = T — strong liniAr^oo S N S N via the 
pulsation uj, see (2.13). It is interesting to observe that the same function ^(S* , S + S~) appears 
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in the computation of L(S + S~) := T — strong lini7v->oo L(S^S^). Again, since (S^S^Y = 
StfSx, we can use formulas (2.32) and (2.33). Here the computations are significantly harder, 
but no difficulty of principle arises. As a technical tool it is convenient to use the fact that, 
in the limit N — > oo, all the intensive operators commute with all the local operators of the 
system, lim^^^S^, <jj ] = 0, for all a, (3 and j. Therefore we get 

L(S + S~) = - 169 ^ S ^ 3 h(S°, S+S-). (3.6) 

The phase structure of the model is now given by the right-hand sides of equations (3.4) 
and (3.6), see [1, 2], and, in particular, from the zeros of the functions 

fi(x,y) = - 89 J 3 V h(x,y), f 2 (x, y) = -^-^-h(x, y), (3.7) 

where we have introduced, to simplify the notation, x = S° and y = S + S~, so that uj = 
gyjx 2 + 4y and v = 2e+gx. In particular, the existence of a super-conducting phase corresponds 
to the existence of a non trivial zero of f\ and / 2 , [1, 2]. Due to the definition of fi and f 2 it 
is clear that any (x Q , y ), with x a ^ and y Q ^ 0, is such that fi(x Q , y ) = /^(^o, y ) = if and 
only if it is a zero of the function h: h(x ,y ) = 0. In order to find such a solution, it is first 
necessary to obtain an explicit expression for the coefficients 9ftr±\ This is easily done using 
the definitions in (2.26), since we get 



£p(a) 



1 r°° 

± -3/ E \UP)\ 2 e- iT1/±(pl dT = n £ \fm(p)\ 2 6(v ± (P)), (3-8) 
and 

*tr$!=* x: \fn(p)\ 2 s(Mp}). (3.9) 

peAjv 

It is now almost straightforward to recover the results of [1, 2]. Following Buffet and Martin's 
original idea, we look for solutions corresponding to v — 0. This means that, because of (2.13), 
the value of x = S° is fixed: x = —2e/g. Moreover, with this choice, v+(p) = uj — ep, which 
is zero if and only if uj = ep. Also, we have v~(p) = —uj — ep, which is never zero. For these 
reasons we deduce that = 0, 7 = a, b, while the sums in (3.8) and (3.9) for ftrj are 

restricted to the smaller set, Sn C A^r, of those values of p such that, if q e £n then e$ = uj. 
Therefore, recalling the expression of m(p) and n(p) in (2.18), we find 

^ = E 3Hf = £ \f(p)\\ (3.10) 
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From definition (3.5), therefore, we get the following equation 

or 

e** = 9 -±^. (3.11) 

This equation is the crucial one, which replaces the equation obtained in [1, 2], g tanh {^f) = w. 
We conclude that: 

1) first of all, introducing a new variable £ = ^, equation (3.11) has a non-trivial solution if 
and only if the function g(£) = — ^| has a zero £ 7^ 0. It is clear that such a solution does 
exist only if the first derivative of <?(£), computed in £ = is positive, i.e. when (3g — 2 > 0. 
This is because g(0) = and lim^x- g(£) = —00. We recover therefore the first result of [1, 2], 
since this inequality implies the existence of a critical temperature, T c := coinciding with 
that found by Martin and Buffet, such that, when T < T c , a £ 7^ does necessarily exist such 
that (?(£) = 0, and the system is in a superconducting phase. 

2) it is also possible to find the value of y = S + S~ directly from equation (3.11). However, 
in order to recover the same value of the gap operator known in literature, we prefer to play a 
little bit with equation (3.11) in the following way: 

(3u\ e~ - e-~ e ft"- 1 ~ 1 

z J ea+ea e^+i — + 1 

This chain of equalities shows once again how our equation (3.11) returns the same equation 
obtained in [1, 2] with completely different techniques. 

Vice versa, it is also straightforward to check that equation g tanh (j^j = oj implies equation 
(3.11): 

tanh(f) + l | + 1 g + LU 
tanh(^)-l 1-f g-u' 
and this concludes the proof of the equivalence of our approach with that of Buffet and Martin. 



IV Conclusions and comments 

We have shown how the SLA can be successfully used to analyze the phase structure of low 
temperature superconductivity analyzing a strong coupling BCS model, considered as an open 
system interacting with a bosonic thermal bath. 
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The procedure, which makes use of the canonical representation of thermal states, is rather 
direct and is technically much simpler than the one used in the original paper, [1] . Among the 
other simplifications, for instance, a single equation h(x, y) — must be solved instead of the 
system fi(x, y) = f2(x, y) = 0, which is the highly transcendental system which must be solved 
in [1]. 

For this reason we believe that it may be worth considering other models, still unsolved, 
with the simplifying strategy provided by the SLA, since new insights may eventually come 
out. For instance, one could first replace the bosonic reservoir with a reservoir made of quons, 
[9], in the attempt of getting a different free time evolution for the creation and annihilation 
quon operators. Following our analysis, and in particular definition (3.5) of the function h, this 
is in fact the easiest way to get an higher value of the critical temperature (T c > ^). Another 
possibility to achieve the same result is to consider a second reservoir interacting with the first 
one: again, in this way the free time evolution of the bosonic operators will be different from 
the one considered here, aj(fe~ tet ). These models will be considered in a forthcoming paper, 
[10]. 
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Appendix: Few results on the stochastic limit 



In this Appendix we will briefly summarize some of the basic facts and properties concerning 
the SLA which are used all throughout the paper. We refer to [3] and references therein for 
more details. 

Given an open system S + 1Z we write its hamiltonian H as the sum of two contributions, 
the free part H and the interaction XHj. Here A is a coupling constant, H contains the free 
evolution of both the system S and the reservoir 7Z, while Hi contains the interaction between 
S and TZ. Working in the interaction picture, we define Hj(t) = e lH o t ff ie - lH ot an( j the so called 
wave operator U\(t) which is the solution of the following differential equation 

d t U x (t) = -iXH^Uxit), (A.l) 

with the initial condition U\(0) = 1. Using the van-Hove rescaling t — > -p-, see [2, 3] for instance, 
we can rewrite the same equation in a form which is more convenient for our perturbative 
approach, that is 

with the same initial condition as before. Its integral counterpart is 

U X (±) = 1 - j f Q Hj^U^dt', (A.3) 

which is the starting point for a perturbative expansion, which works in the following way. 

Suppose, to begin with, that we are interested to the zero temperature situation. Then let 
y?o be the ground vector of the reservoir and £ a generic vector of the system. Now we put 
ip^ — tp <S> £. We want to compute the limit, for A going to 0, of the first non trivial order of 
the mean value of the perturbative expansion of U\(t/\ 2 ) above in tpffi, that is the limit of 

h(t) = {-\) 2 Jl dt^dh <H I ( t ^)H I ( t ^) >^), (A.4) 

for A — > 0. Under some regularity conditions on the functions which are used to smear out the 
(typically) bosonic fields of the reservoir, this limit is shown to exist for many relevant physical 
models, see [3], and [7, 6] for few recent applications to quantum many body theory. It is at 
this stage that all the complex quantities like the T^'s we have introduced in the main body 
of this paper appear. We define I(t) = lim^o I\ (t) ■ In the same sense of the convergence of 
the (rescaled) wave operator U\(-k) (the convergence in the sense of correlators), it is possible 
to check that also the (rescaled) reservoir operators converge and define new operators which 
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do not satisfy canonical commutation relations but a modified version of these. For instance, 
in Section II this procedure has produced the operators c^ 7 ] starting from cpj. Moreover, these 
limiting operators depend explicitly on time and they live in a Hilbert space which is different 
from the original one. In particular, they annihilate a vacuum vector, 770, which is no longer 
the original one, tp . This is what happens, for instance, if tp depends on A, </? — > V?o*\ and 
considering 770 as the following limit: 770 = liniA^o V 9 ^ • 

It is not difficult to deduce the form of a time dependent self-adjoint operator Hj Sl \t), which 
depends on the system operators and on the limiting operators of the reservoir, such that the 
first non trivial order of the mean value of the expansion of Ut — 1 — % Jq Hj Sl \t')Ut'dt' on the 
state 77^ = r] <g> £ coincides with I(t). The operator U t defined by this integral equation is 
called again the wave operator. 

The form of the generator follows now from an operation of normal ordering. More in 
details, we start defining the flux of an observable X = X <g> l r , where l r is the identity of the 
reservoir and X is an observable of the system, as jt(X) = U}XU t . Then, using the equation 
of motion for U t and [//, we find that d t j t (X) = iU}[H\ sl \t), X]U t . In order to compute the 
mean value of this equation on the state tjq , so to get rid of the reservoir operators, it is 
convenient to compute first the commutation relations between Ut and the limiting operators 
of the reservoir. At this stage the so called time consecutive principle is used in a very heavy 
way to simplify the computation. This principle, which has been checked for many classes of 
physical models, [3], states that, if (5{t) is any of these limiting operators of the reservoir, then 

\P(t),U v \ = 0, for all t > t'. (A.5) 

Using this principle and recalling that rjo is annihilated by the limiting annihilation operators 
of the reservoir, it is now a simple exercise to compute < dtjt(X) > to and, by means of the 

Vq 

equation < d t j t (X) > (o=< jt{L{X)) > (o , to identify the form of the generator of the physical 
system. 

Let us now consider the case in which T > 0. In this case the state of the reservoir is no 
longer given by ip . It is now convenient to use the so-called canonical representation of thermal 
states, [3]. Using the same notation of Section 2, any annihilator operator be written 

as the following linear combination 

or,- = V^cg + yfi(0j<$\ (A.6) 
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where m{p) and n{p) are the following two-points functions, 

. 1 , e -/3ep 

m(p) = Ma^alJ = ra(p) = ^(ata fo -) = j-— (A.7) 

for our bosonic reservoir, if up is a KMS state corresponding to an inverse temperature (3. The 
operators are assumed to satisfy the following commutation rules 

[cg.cgVW, (A.8) 
while all the other commutators are trivial. Let moreover $ be the vacuum of the operators 

(a) 
C - ■ '. 

cg$ = 0, Vp,j,a. 

Then it is immediate to check that the results in (A.7) for the KMS state can be found, using 
these new variables, representing uop as the following vector state u>p(-) =< $0) '^o >• With this 
GNS-like representation it is trivial to check that both the CCR and the two-point functions 
are easily recovered. This representation is also called in [3] the Fock-anti Fock representation 
because of the different sign in the free time evolution of the annihilation operators c^j and 
Cpj. Once this representation is introduced, all the same steps as for the situation with T = 
can be repeated, and the expression for the generator can be deduced using exactly the same 
strategy. 
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